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A generalization of the Van der Waals excluded volume procedure for the multicomponent hadron 
gas is proposed. The derivation is based on the grand canonical partition function for the system 
of particles of several species interacting by hard core potentials. The obtained formulae for ther- 
modynamical quantities are consistent with underlying principles of statistical mechanics as well as 
with thermodynamical identities. The model can be applied to the analysis of experimental data 
for particle number ratios in relativistic nucleus-nucleus collisions. 
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^ \ I. INTRODUCTION 

Thermal hadron gas (HG) models have recently been used to fit the data of particle yields in nucleus- nucleus 
(A+A) collisions at the AGS and SPS energies (see e.g. [Q). The ideal HG model becomes inadequate at the chemical 
Cn , freeze-out in high-energy A+A collisions: temperature and baryonic chemical potential obtained from fitting the 
particle number ratios at the AGS and SPS energies lead to artificially large values of total particle number and 
energy densities (see e.g. Q). This is hardly consistent with a picture of a gas of point-like noninteracting hadrons. 

The Van der Waals (VdW) excluded volume procedure appeared to be effective in taking into account the hadron 
\Q . repulsion at short distances. It suppresses undesirable large values of particle number densities. Different versions 
of the VdW HG models were proposed and applied for fitting experimental data on particle number ratios in A+A 
collisions at the AGS and SPS energies [3-12]. The proper volume of the i-th hadron species is expressed in terms of 
its hard-core radius Ri. An introduction of the phcnomcnological parameters Ri changes the particle number ratios 
in comparison with ideal HG results. The VDW model formulation, however, has not been properly defined in the 
case when the RiS are not equal to each other. 

The aim of the present paper is to propose a generalization of the VdW excluded volume procedure for the HG gas of 
several particle species with different particle radii. The derivation is based on the grand canonical partition function 
for the system of particles of several species interacting by hard-core potentials. The obtained formulae are therefore 
consistent with the underlying principles of statistical mechanics as well as with thermodynamical identities. The 
pressure, particle densities and other thermodynamical quantities as functions of temperature and chemical potentials 
are defined by a set of coupled transcendental equations. 
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II. ONE COMPONENT VDW GAS 

The canonical partition function (CPF) for the one-component classical (Boltzmann) gas can be written as 

where V and T are the system volume and temperature, m and N are the mass and number of particles, respectively. 
The function U in Eq. (|lj) is assumed to be equal to the sum of pair potentials: 

U = ]T «(|fi-fi-|). (2) 

l<i<j<N 



After integration over the particle momenta, Eq. (TO) is reduced to 
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Here wc use the notation 
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where Ki is the modified Bessel function. The asymptotic behavior of <fi(T;m) in the non-relativistic limit, m >> T, 
has the form 



<p(T;m) 
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exp(— m/T) . 



By means of the Mayer functions 



exp 



T 



one can rewrite the integrand of Eq. (|3j) in the following form 

, TJ \ N—l N jv-i iv 
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The approximate equality in the last expressions is based on the assumption that the gas is rarefied and, therefore, 
the terms containing products fijfim can be dropped. 
Let us introduce the notations 



d 3 nf %3 = - 2b{T,?j) 



Now we shall use the rigid ball model, i.e. we assume the hard-core interaction between the particles: 

udn-fjl) : 



oo if \fi -fj\ < 2R 
if \Fi-rjl >2R, 
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(R is the particle radius). In this case b defined by Eq. (0) does not depend on temperature. For T/ 1 / 3 ;§> R the 
dependence of b on fj is negligible and the integral in Eq. (pi) can be easily calculated: 



16 , 

b = y ^3 



Substituting Eqs. (0]|) into Eq. (||) one gets 

1 N 

Z(V,T,N) = [ (f) (T;m)} N l[(V-2(N-i)b) 
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In a rarefied gas the value of 2bN is much smaller then the total volume V. In this case one can approximate 
log(l -x)^-x in Eq. (pi). This yields 



Z(V,T,N) 
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Performing the summation in the exponential and using approximate relation exp(— Nb/V) = 1 — Nb/V one gets 
CPF of one component rigid ball gas: 



Z(V,T,N) - jL[4>(T;m) V\ N exp(-N(N-l)^ 
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Substituting the expression ( |13[ ) into the formula for pressure one finds the well known VdW equation for rigid ball 
gas: 

„ )sT ^^>=_^_. (14) 

Note that b ( |l0| ) is equal to Av, where v — 4nR 3 /3 is the particle volume. The VdW equation of state (|14| ) is obtained 
from the statistical mechanics of rigid balls within a rarefied gas approximation, vN << V. For a dense gas the VDW 
equation (y_4|) should be considered as a phenomenological extrapolation. 

The grand canonical partition function (GCPF) is expressed through the CPF in the following way: 

Z(V, T, M ) = ]T exp (£p) Z(V, T, N) , (15) 

N=0 ^ ' 

where /i is the chemical potential. In the case of rigid ball gas, the upper limit of the sum is in fact not infinite. The 
CPF (pi) becomes equal to zero if N exceeds Nq ~ V/v. 

Let exp I LLjr- J Z(V,T, N*) is the maximal term in the sum (15). The rest of terms are also positive, hence the 
following inequalities are satisfied 

e*p(^) Z(V,T,N*) < Z(V,T,n) < N x exp (^\ Z(V,T,N*) , (16) 

The pressure can be expressed via GCPF by the formula 

P{T,fi) = T hm . (17) 

v— >oo v 

Taking Eq. (Ru) into account one gets 

(T , . r yN*+T log Z(V,T,N*) 
p(T,n) > hm , 

(T x . r fiN*+T log Z(V,T,N*) logJV 

p(T,n) < hm + T hm — — - . 

Since V^ 1 log N — > in the thermodynamical limit, V — * oo, the pressure p(T, /j,) ( p7| ) is defined by the largest term, 
with N — N*, of the GCPF ( jig). N* is also the average number of particles in the grand canonical formulation. 
Using the VdW approximation (|l3|) for the CPF one finds 



P{T,n) = T lim — log 

V— »oo V 
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where A — exp(/i/T)0(T; m), and iV* = N*(V,T, fi) corresponds to the maximum of the expression in the square 
brackets. 

Let us show that Eq. ( |l8| ) leads to the result 

p(T,n) = T£, (19) 

where £ is defined by the transcendental equation 

£ = A exp(-&0 . (20) 

Using the asymptotic representation for the T-function logarithm at N — * oo 

logr(JV+l)S JV(logiV-l) (21) 

it is easy to check that the value of N* satisfying the maximum condition of logarithm argument in Eq. (ITq) is given 
by the formula 

N* =* Vn , (22) 



where n = n(T, /i) is related to £ via equation 



l + 6£ 



(23) 



Substitution of Eq. ( P2|) into Eq. ( ^8| ) with account for Eqs.(2fj2J|) yields the formula (|l9|). The quantity n = n(T, //) 
is a particle number density in the grand canonical formulation. One can readily check that the definition of the 
particle number density, n = dp(T, fJ.)/d/j, , leads to Eq. (p3|), provided that Eqs. (|l9|-f2l|) are taken into account. 
For the point like particles, R = and 6 = 0, Eq. ( |19| ) is reduced to the ideal gas result: 

p id (T,^) = T exp(/i/T) <j>(T;m) = T n ld (T^) . (24) 

Eq. ( |l9| ) can be therefore written in the form j|: 

p(T,,)=^(t,,-^-) . (25) 

It can be also presented in the form of Eq. ( jl4| ) with N — N*(V, T, /i), which demonstrates explicitly the equivalence 
between canonical and grand canonical formulations at V — > oo. 

III. TWO COMPONENT VDW GAS 

In the case of two particle species CPF has the following form: 



X II (0^3 GXP ~ T eXP 



(2tt) 3 M r V r 

fc=i > ' 

where mi, N± (m2, iVjj) are the mass and number of particles of the 1-st (2-nd) species, 

U™ = £ «n(|rf ) -rf|)+ £ ^(Irf-rfl) (27) 

l<i<i<ATi l<k<l<N 2 

iVi JV 2 

1=1 fc=l 



vyv~< ^" - " 2i 



It should be mentioned that in the case of two particle species U^ 1 ' 2 ^ contains three types of two-particle potentials. 
While the potentials itn and U22 describe interactions between particles of the same species and can be handled 
similarly to the potential u of one-component case, the potential u\2 describing interactions between particles of 
different species requires an special treatment and prevents a straightforward generalization of the one-component 
VDW equation to the two-component gas. 

The integration of ( |26| ) over the particle momenta gives the following expression for the CPF 

[0(T; mi) r \m m 2 )} N * f g d 3 r (D -Q ^.CO _ f ^ ' 

i=i fe=i 

The next step, however, involves the Mayer functions of three types (p, g=l,l;2,2;l,2): 



ww ^^ff^ / n^fl'-M-^F) ■ < 28 » 



/'r' = expl- """'- r '' " I - 1, p., = 1,2 (29) 

(note that f} 2 — /y 1 ). The integrand of (ES) can be rewritten as 
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where u, w = 1,2 and iV w > N w . 

For each type of Mayer function we introduce the notation 



JdM p) f^ = -2b pg (T^). 



(31) 



For the rigid ball model at V 1 ^ 3 ^> max{R\, R 2 ) (R P is the radius of a particle of p-th species, p = 1,2) it yields 

2 
b Pi = T^{R P + RgY 

Substituting @ and @ into @ one gets 

Z(V,T, #!,#;,) = j^j^ [ct>(T;mi)} Nl [<f>(T;m 2 )] N2 f[(V - 2(N W - k)b ww ) 
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We again assume that the gas is rarefied and the total proper volume of all particles is much smaller then the total 
volume of the system. Then it follows from Eq. (p3|) : 



Z(V,T,N X ,N 2 ) 
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Performing the summation in the exponential 

1 



Z{V,T,N U N 2 ) 
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and introducing quantities b\ 2 and 6 2 i constrained by the condition 

b\ 2 + hi = 2&12 , 
we can rewrite the CPF in the form 

Z{V, T, N lt N 2 ) 2 jJ^j [4>{T- m,)]" 1 [<f>{T; m 2 )f* V N ^+ N ^ 



(35) 



(36) 



exp I -Ni— ~ N 2 — 
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exp I -N 2 — -iVi — 



A r 2 



(37) 



Imposing the additional constraints 2N v b vw /V « 1 we obtain the final expression for the partition function of the 
two-component VdW gas 



Z(V,T,Nx,N 2 ) 



1 



[4>(T; mi )] Nl [<KT;m»)f> 
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Substituting expression (pq ) into the formula for the pressure one obtains 

,siog^(v;r,JVi,JV2) n ± t 



p(V,T,N u N 2 )=T- 



N 2 T 



dV 



V - b 11 N 1 - b 21 N 2 V - b 22 N 2 - b 12 Ni 



(38) 



(39) 



The equation of state (B9J) is consistent with the virial expansion up to second order and can be considered as a 
generalization of the VdW equation for the two-component system. We define 612 and 621 as the following: 



615 



61161 



611 



621 



62261 



611 + 6 2 



(40) 



Eq. (J40j) satisfies the constraints (|36|) and leads to the correct physical behavior in the limiting cases when both particle 
radii are equal to each other or when one of the particle radius is equal to zero. If two species have equal radii, R\ = R 2l 
equation ( |39| ) is reduced to the one-component VdW equation ( fL4[ ) with N = Nx + N 2 and 6 = 611 = 622 = 612 = 621. 
Note that in a general case, R\ ^ R 2 , the VdW excluded volumes in Eq. (J3^) are different for different particle 
species. 

The transformation to the grand canonical ensemble is similar to that of one-component case. The GCPF has the 
form 



00 00 



Z(V,T,fxi,n 2 ) 
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(41) 



where fi q (q = 1, 2) are the chemical potentials of each particle species. The total number of terms, No, in double sum 
( [4l| ) is of the order of V 2 . Since V^ 1 log No — > at V — + 00 the pressure can be expressed via the maximum term in 
the sum @): 



p(T,m,[Mi) 



T 
lim — log 

y^oo V 



exp 
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In the VDW approximation (|38|) the last expression takes the form 

,ac a n: 
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p(T,m,n 2 ) = lim — log 
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where A p = cxp(/j p /T) 4>{T,m p ) . 

Let us show that the pressure (B3h can be calculated by the formula 

P(T, Hi, H2) =T(£i + £ 2 ) , 
where the values of £ g are found from the set of coupled transcendental equations 

£1 = A x cxp(-6n£i - 6126) > 
£2 = M exp(-6 2 26 - 6 2 i£i) . 



(44) 



(45) 
(46) 



Using the asymptotic representation for the T-function logarithm it is easy to check that the values of N* satisfying 
the maximum condition of the logarithm argument in Eq. (J43j) are given by the formula 



N* p 2 Vn p , 

where n p = n p (T,fXx,jJ< 2 ) are related to £ p via the equations 

nx 



£1 = 

£2 = 



1 - nxbxx - n 2 b 21 
H2 

1 - n 2 b 22 - nxb 12 



(47) 

(48) 
(49) 



The substitution of Eq. (E7|) into Eq. (E3|) yields the formula (E 



The set of linear equations (48) and (}49j) can be solved for n p : 



m = = — = — , (50) 

1 + fl&ll + 6^22 + 66(6ll6 2 2 - 612621) 

6[i + 6(611 -612)] r ,.s 

n 2 = = — = . (51) 

1 + £1611 + £2622 + 66(611622 - 612621) 

Similarly to the one component case, the quantities n p are particle number densities. One can readily check that the 
definition n p — dp(T, fix, fi<i)/dfi p leads to Eqs. ( |50| ) and (pl|), provided that formulae (^446[) are taken into account. 

IV. MULTICOMPONENT VDW GAS 

The above considerations can be generalized to the multi-component VdW gas with an arbitrary number of particle 
species. After integration over particle momenta and simplifications similar to those of the above consideration one 
gets the expression for the CPF of if-component VdW gas: 

k / k \ N « 

Z(V,T,N 1 ,...,N K ) ^ ]J — [cb(T; mq )}^[v~Y,N P b pq \ , (52) 

where 

- _ 2b pp b pq 
pq ~ b +b ' [ ' 

PP ^ 1Q 

2 
b pq = b qp = -n{Rp + R q ) 3 . (54) 

R q is the radius of particle of speciesg. 4>(T;m q ) is defined by Eq. (Eh. The CPF ( p2[ ) yields the VDW equation of 
state for a A"-component rigid ball gas 

K NT 

P{V,T,N X ,...,N K ) = £ vX ,k l N - (55) 

q=\ V l^p=l °pq ly p 

The pressure in the grand canonical ensemble is given by the formula 

K 

p(T, [iu -, m) = T Y^p ' (56) 

P =i 

where \i v are the chemical potentials (p = 1, ..., K). The functions £, q satisfy the set of coupled transcendental equations 

(57) 



6 = A p exp I - Y kpqtq ) , 



where A p = exp(fi p /T)(j)(T;m p ). The particle densities n p = n p (T, fix, ...,/iR-) are obtained as the solutions of the 
following set of coupled linear equations 

( > - i-jf"„i ■ (58) 

1 L-iq=\ ,i q u qp 

In the canonical ensemble formulation of the HG model the numbers N±, ..., Njc are not fixed. They can not be fixed 
because of inelastic reactions between the hadrons. The fixed values have the conserved charges: baryonic number 
B, strangeness 5 (strangeness is conserved as we neglect weak decays) and electric charge Q. The CPF has then the 
form 



oo K 

Z{V,T,B,S,Q) = Y, Z(V,T,N X ,...,N K ) 5{B -J^W (59) 

N 1 ,...,N K = l t=l 

K K 

x *(S-5^s«JVi)5(Q-5^©JVi), 

t=l i=l 

where 6j, s, and qi are the baryonic number, strangeness and electric charge of i-th hadron species. In the application 
of the thermal HG models to A+A collisions, B equals the number of nuclcons participating in the reaction, 5 = 
and Q — aB with a « 0.5 for intermediate nuclei and a ~ 0.4 for heavy nuclei. CPF ( |59| ) can be calculated with the 
VdW input (|5^) for Z(V, T, N\, ..., Nk)- The presence of the Kronecker (5-functions in Eq. (|59|) makes the canonical 
ensemble rather complicated. The grand canonical formulation justified for large systems is much more convenient. In 
this case the system properties are defined by the pressure function ( pq ) with the chemical potentials fit (i = 1, ..., K) 
defined as 

Hi = bijX B + SifXS + Qi^Q (60) 

in terms of the baryonic fiB , strange us and electric ^iq chemical potentials. They are chosen, at given V and T, 
to fix the average values B, S = and Q = aB. 

V. SUMMARY 

In the present paper we have proposed a generalization of the VdW excluded volume procedure for the multicom- 
ponent hadron gas. The canonical and grand canonical ensemble formulations are presented. Different hard-core radii 
for different hadron species have been discussed in the literature (e.g. [8-12]). The excluded volume procedure was 
based on the substitution V — > V — J2i bu^i- This ansatz is shown to be not correct: for unequal hard-core radii the 
VdW excluded volumes seen by the different hadron species are different. Thermodynamical quantities in the grand 
canonical formulations are found to be the solution of a coupled set of the transcendental equations. 
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